Pair-tunneling resonance in the single-electron transport regime 
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We predict a new electron pair-tunneling (PT) resonance in non-linear transport through quan- 
tum dots with positive charging energies exceeding the broadening due to thermal and quantum 
fluctuations. The PT resonance shows up in the single- electron transport (SET) regime as a peak 
in the derivative of the non-linear conductance, d?I /dV^ , when the electrochemical potential of one 
electrode matches the average of two subsequent charge addition energies. For a single level quan- 
tum dot (Anderson model) we find the analytic peak shape and the dependence on temperature, 
magnetic field and junction asymmetry and compare with the inelastic cotunneling peak which is of 
the same order of magnitude. In experimental transport data the PT resonance may be mistaken for 
a weak SET resonance judging only by the voltage dependence of its position. Our results provide 
essential clues to avoid such erroneous interpretation of transport spectroscopy data. 

PACS numbers: 73.63.Kv, 85.35.Gv, 85.35.-p, 
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Introduction. Electron tunneling spectroscopy has 
nowadays become a standard tool for investigating the 
in situ properties of single-electron transistors based 
on quantum dots in semi-conductor hetero-structures, 
nanowires, carbon nanotubes and even single molecules. 
The basic spectroscopy rules derive from the simple en- 
ergy resonance conditions for single-electron tunneling 
(SET) onto the dot. As a result, the bias positions of 
differential conductance (dl/dV) resonances depend lin- 
early on the gate voltage due to capacitive effects [l[ . In 
addition, inelastic cotunneling (COT) processes can of- 
ten be distinguished as steps in dl/dV 0, allowing for 
spectroscopy of quantum dot excitations with increased 
energy resolution Q . Here an electron is transferred from 
the high to the low biased electrode through the dot, us- 
ing the excess bias-energy to excite the dot by an energy 
Ae. Since these electron-hole charge transfer resonances 
involve only a virtual charging of the dot, they appear at 
a bias threshold V = Ae, independent of the gate voltage. 
These processes arise only in second order perturbation 
theory in the tunnel rate F. 

In this letter we show that, in the same order of F, elec- 
tron pair-tunneling (PT) gives rise to distinct, measur- 
able transport effects which, to our knowledge, have been 
overlooked so far. In a single coherent process, a pair 
of electrons is extracted from (added to) one electrode, 
resulting in a real (dis)charging of the dot by two elec- 
trons, involving a positive charging energy much larger 
than broadening due to thermal and quantum charge 
fluctuations. The corresponding PT resonance appears 
deep inside the SET region where a finite current is flow- 
ing and its position has the same gate-dependence as 



SET resonances. It might thus be mistaken for a weak 
SET resonance belonging to some excited state and one 
might thereby extract an erroneous level-structure from 
the spectroscopic data. We show how the PT resonance 
can be distinguished from other types of first and second 
order processes and what additional information can be 
extracted from it. Since spectroscopic data obtained from 
single-electron transistor devices are often non-trivial to 
interpret we provide several independent criteria for its 
experimental identification, namely the resonance shape 
and the dependence on temperature, magnetic field and 
junction asymmetry. The PT processes discussed here 
completely lack any signature inside the Coulomb block- 
ade region where SET is exponentially suppressed. In 
this region a different type of pair-tunneling resonance 
occurs which is much weaker than the one discussed here 
since it appears only in the third order in the tunnel 
rate F Also the latter pair-tunneling process in- 

volves neither real nor virtually doubly occupied states 
(i.e. remains finite in the U = oo limit). Second or- 
der pair-tunneling resonances were discussed previously 
for the rare case of negative Coulomb charging energy 
(i.e. effective attractive electron-electron interaction) 5|, 
related to pair-tunneling in superconducting grains 6|, 
where SET is suppressed by the superconducting gap. 
Also Kondo effect involving electron pairs have been dis- 
cussed in this limit . Here we consider the experimen- 
tally most frequently occurring case of strong positive 
charging energy. 

Model. The pair-tunneling resonance is already present 
in the generic model for a quantum dot, i.e. the non- 
equilibrium Anderson model. Here the quantum dot. 
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FIG. 1: (Color online). Energy differences between many- 
body eigenstates (dot chemical potentials) of the Anderson 
model and sketch of a COT process (a) and a SET process 
(b). (c): More detailed sketch of a PT process, giving rise to a 
resonance at A = [7/2 above the SET onset. Filled (unfilled) 
circles indicate real (virtual) occupation. 



described by the Hamiltonian H = eo-ncr + Ufi-^hi, 
consists of a single orbital with energy and occupa- 
tion fia = d\da- Here e| — e| = h equals the Zee- 
man energy and U is the finite positive charging en- 
ergy. We let N = her denote the electron number 
on the dot. The many-body eigenstates are |0), and 
2) with energies = 0, E'er = e^- and E2 ~ + U. 

The dot is coupled by a tunneling Hamiltonian. ~ 
J2r=L,RJ2k.a'^rdlcrka- + h.c, to macroscopically large 

reservoirs, described by Hb. = J2r=L.R^k,cy ^kcli^^Crka- 
The electrons in the reservoirs are assumed to be non- 
interacting, with operators c^^j.^., Crka for state k and spin 
a in electrode r = L,R. The tunnel amplitudes can be 
expressed in the tunnel rates Tr through Tr = y^Tr/iTT. 
Throughout the paper we use natural units where h — 
= |e| = 1, where — |e| is the electron charge. 
Let us first present the basic physics. Fig. [T] shows 
a schematic comparison of tunneling processes for the 
spin-degenerate Anderson model when E^r is larger than 
the average chemical potential n of the reservoirs, such 
that the quantum dot is unoccupied at zero bias voltage, 
V = 0. We take fi = and assume symmetric biasing 
such that /XL = 1^/2, /^R, = — y/2. At low bias transport 
is dominated by cotunneling involving virtual occupation 
of state |cr). For larger bias, \V/2\ > E^ — Eq, electrons 
can sequentially tunnel into and out of the dot, involving 
real occupation of state \a). In this regime double occu- 
pation of the dot through two consecutive SET processes 
only becomes energetically allowed when additionally the 
energy difference E2 — E^ is below the chemical potential 
of one lead, i.e. when \V/2\ > E^ -\- U . Midway between 
these resonances the coherent tunneling of a pair of elec- 
trons from the same reservoir becomes possible, i.e. at 
the resonance condition: 



\V/2\>\{Y.„E„ + U). 



(1) 



The corresponding process is shown in Fig.[T]Jc). One can 
think of this as one electron from just below the Fermi 
level of the left reservoir tunneling onto state jcr), leaving 



an excess energy A = C//2, which can be used to assist 
the second electron in reaching state |2). Thus, the total 
PT process is energy-conserving. However, since there 
are no internal degrees of freedom on the dot to store the 
energy A, these two processes have to take place coher- 
ently in the short time set by the time-energy uncertainty 
relation. 

Experimentally, the dot energies can be controlled lin- 
early by the gate voltage Vg. This produces the charac- 
teristic Coulomb diamond figures in differential conduc- 
tance color- maps as function of gate- and bias voltage [l[ . 
Since the resonance condition ([T]) for PT depends on the 
average of the two dot energies E„ the PT lines have 
the same slope as SET resonances. To distinguish the 
PT and SET resonances a consistent calculation of their 
transport signature is thus imperative. 

Transport calculation. The pair-tunneling resonance 
is visible in the transport spectrum when the charging 
energy is larger than both the tunnel rate F and temper- 
ature T, and originates from coherent two-electron pro- 
cesses which become important with increasing F. We 
can address the interesting regime U > T > T using the 
real-time transport theory which allows systematic 
treatment of processes beyond lowest order in F, while ac- 
counting for the competition of single- and two-electron 
processes, essential at the PT resonance. The central 
task is to calculate the stationary reduced density ma- 
trix of the dot in the basis of many-body eigenstates, 
which is done using a kinetic (master) equation. From 
this density matrix the current flowing out of reservoir r, 
/r, can be calculated. The required transport rates are 
calculated perturbatively up to second order in F and are 
given by the analytic expressions which we derived for a 
general model in Ref. [9|. Importantly, all coherent one- 
and two-electron processes are included and the Coulomb 
interaction U is treated non-perturbatively. 

Fig. [2] (a) and (b) show the differential conductance as 
a function of gate and bias voltage, where in (b) the spin- 
degeneracy is lifted by an applied magnetic field. SET 
processes give rise to the strong yellow lines which cross 
at zero bias at the charge degeneracy points. In (b) ad- 
ditional gate-independent steps show up in the = 1 
Coulomb blockade region, indicating the onset of inelas- 
tic cotunneling leading to occupation of the excited spin- 
state. Subsequent relaxation of this state by SET gives 
rise to additional gate-dependent lines in the Coulomb 
blockade region, as described e.g. in Refs. 111- 
Clearly separated from these other resonances, the PT 
resonance associated with an electron pair tunneling onto 
(off) the dot appears as a gate-dependent step inside the 
SET region, indicated by a red sohd (green dotted) ar- 
row in (a). Considering the resonance position and its 
voltage dependence only, the PT could be mistaken for a 
weak SET process involving an excited state of the = 1 
charge state with energy A = U/2 relative to the A^ = 1 
ground state. However, the resonance shape provides one 
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FIG. 2: (Color online), (a): dl/dV vs V,Vg plotted on log- 
arithmic color scale for Tl = Tr = T, = lOOT = 2000r 
and Cr = Cl = lOOCg. (b): Same as (a) except for the 
finite Zeeman splitting hj h = 15T and asymmetric capac- 
itances Cl = 0.7Cr = lOOCg. In addition, the gate capac- 
itance of state 2 is larger: Cg,2 = 1.176Cg. (c): d?I/dV^ 
vs V around the PT threshold Vpt at Vg = -0.625C/C/Cg 
(C = Cl -I- Cr + Cg), indicated by the red line in (a). This 
demonstrates the scaling with l/U"^ and the agreement with 
the analytic expression Eq. (|3} marked by (A), (d): The 
two independent Keldysh diagrams contributing to the pair- 
tunneling resonance. A sum over a is assumed and a denotes 
the opposite spin-projection of cr. 



crucial clue to the identification as PT: in Fig. [5Jc) we 
show a bias trace along the red vertical line indicated in 
FigHUa), where it is seen that the PT resonance appears 
as a peak in cfl/dV^, rather than in dl/dV as for SET. 
The reason is that, analogous to inelastic cotunneling, 
the number of electron states in the reservoirs available 
for a PT process is proportional to the bias voltage [3]. 
Another clear way to separate PT and SET resonances 
is that the PT resonance exhibits no Zeeman splitting, 
see Fig. ^h) . The physical reason is that the tunneling 
of a pair involves a transition between states with the 
same spin (iV = 0,2 electron singlets), expressed by the 
appearance of the sum E-^ + Ei in Eq. ([1]). Since PT 
relies on an empty or doubly occupied initial state, the 
resonances do not continue down into the linear trans- 
port regime, making PT an exclusively non- equilibrium 
effect. Additionally, in contrast to the magnetic field ex- 
citations, the PT resonances do not continue horizontally 
into the N ~ 1 Coulomb blockade region as inelastic co- 
tunneling steps. In cases which deviate from the simplest 
capacitivc model, the PT resonance shows an additional 
distinct property. In Fig. [DJb) the capacitances associ- 
ated with the left and right leads were chosen unequal 
and the doubly occupied state has a larger gate-coupling 



(this might happen if the many-body wave- function of 
the doubly occupied state is localized closer to the gate 
electrode, or is less screened by the source and drain elec- 
trodes [l^)- This causes the diamond to be both tilted 
and skewed, but the PT resonance can still be found by 
taking the average bias voltage positions (c.f. ((1])) of the 
N = ^ N = 1 and = 1 ^ = 2 ground state 
transitions respectively. 

The precise amplitude and shape of the PT peak are 
also crucial for its correct identification. Here an analytic 
expression for the peak is helpful, which can be obtained 
for [/ > T, F, /i. For this we focus on the upper left SET 
region in Fig. [2Ja-b) {V > 0) where the PT involves the 
transition |0) |2) with both electrons tunneling onto 
the dot from the left reservoir. Deep in the SET region 
the SET rates are constant, while all second order rates 
except the pair-tunneling are approximately linear in bias 
voltage. Here the PT peak is given by 
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There are only two independent terms in the perturba- 
tion expansion which contribute to the non-linear voltage 
dependence of the rate for tunneling with an electron pair 
from the left reservoir, W^i^ = Rc{X^^ + -Djp )• These 
terms correspond to the real-time diagrams [8| shown in 
FigElJd). The corresponding integrals can be solved an- 
alytically without the need for ad hoc regularization. 
Eq. ^ shows that the peak in the derivative of the dif- 
ferential conductance directly maps out the energy de- 
pendence of the pair-tunneling rate. Neglecting all terms 
not contributing to the PT peak in the second deriva- 
tive, assuming equal capacitances to the left and right 
reservoirs, we are left with 
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The resonance condition ([T]) is seen in the argu- 
ment of the function F{x) = d^ {xb{x)) /dx"^ = 
(f coth %-\) I sinh^ (f ) , with b{x) = 1/ (e^ - 1) being 
the Bose function and F{Q) — 1/6. 

It is of interest to compare the PT with the inelas- 
tic COT, which is theoretically and experimentally well- 
studied. In general the PT and COT peaks are of com- 
parable magnitude and their dependence on T and U 
are the same, implying that the PT and COT should be 
simultaneously experimentally accessible. A first differ- 
ence is that the PT peak is completely insensitive to re- 
laxation processes, since it involves a transition between 
ground states in different charge sectors. In contrast, the 
COT resonance consists of two parts. The part which 
remains even in the limit where the voltage dependence 
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of the occupations is negligible (relaxation faster than 
the electron cotunneling and / or asymmetric junction 
^L,R » Tr^l), is given by: 



F 



h-V 
T 



(4) 



This expression is valid at the particle hole symmetric 
point (center of the = 1 Coulomb blockade regime) in 
the limit U :^ h > T > T and ^ > 0. For slow relax- 
ation and Fl ~ Fr the COT peak-height is additionally 
increased and also an asymmetry between V > h and 
V < h appears due to the voltage dependence of the oc- 
cupations [l^. For very small fields and Fl — Fr, the 
COT peak height is seen to be given by 9/2 times that 
of the PT peak. A second difference is that the FT peak 
is strongly suppressed for very asymmetric tunnel rates, 
proportional to the square of the smaller coupling and, 
characteristically, the forward and reverse bias PT peak- 
heights differ by a factor 4, in contrast to the COT peaks 
which remain symmetric. Finally, although a magnetic 
field does not cause a Zeeman splitting of the PT peak, 
Eq. ([3]), it does increase the peak height, in contrast to 
the COT peak, Eq. where the situation is reversed: 
the peak height is independent of the field, while the po- 
sition {V = h, independent of V^) is shifted linearly. 

Quantum dots with comparable level spacing and 
charging energy are especially well suited for observa- 
tion of PT since the background SET current is feature- 
less. In the case of a dense excitation spectrum, e.g. 
due to coupling to a localized vibrational mode, addi- 
tional PT resonances complicate a high-resolution spec- 
troscopic analysis. Extending the calculations of the non- 
equilibrium Anderson-Holstein model in Ref. to finite 
U we find that PT resonances give a signature similar 
to SET associated with a very weakly coupled second 
vibrational mode. Thus, the results presented here are 
of general importance for the analysis of models deal- 
ing with the various complexities of realistic quantum 
dots. In the non-interacting limit, U = 0, the PT is 
completely suppressed. In general, PT involves inter- 
ference between the two contributions (diagrams) shown 
in Fig. Eljd). This interference may be constructive, as 
in the Anderson model, or destructive in other models. 
Interestingly the relative sign of the interfering contri- 
butions is sensitive to the spin of the involved quantum 
dot states. For a generalized quantum dot model where 
the doubly occupied ground state is a triplet, we find 
that the PT is suppressed due to quantum interference 
and can be induced by a magnetic field. Finally, the 
PT also shows up in more complex transport quantities 
such as the current noise, which is sensitive to the ef- 
fective charge transferred in the tunnel process (which is 
±2e here). Extending the real-time approach to the cal- 
culation of shot-noise 14, 1^, including non-Mar kovian 
effects, we find that the PT resonance is associated with 
an increased Fano factor. The above examples serve to 



illustrate that pair-tunneling effects are of general impor- 
tance for transport through nanosystems. 

Discussion and conclusion. In this letter we have the- 
oretically predicted a signature of coherent tunneling of 
electron pairs in the single-electron transport regime, for 
quantum dots with positive charging energies exceeding 
the broadening by thermal and quantum fluctuations of 
charge. Current low-temperature measurements can ac- 
cess such fine details in the first three derivatives of the 
current with respect to voltage without dropping below 



the noise level [12| . Also, many molecular and carbon 



nanotube devices couple strongly to the electrodes, mak- 
ing first and second order transport features compara- 
ble. We mention the possibility of measuring the charg- 
ing energy by identifying the PT resonance and reading 
U/2 off in the stability diagram, thus limiting the volt- 
ages one needs to apply. Alternatively, the pair tunneling 
resonance provides a consistency check on SET level as- 
signments. On a general level, the results indicate the 
importance of complete perturbative treatment of non- 
equilibrium problems Q • 
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